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Abstract 



O 

<D ' Impact processes of nanoclusters subject to thermal fluctuations are inves- 

tigated, theoretically. In the former half of the paper, we discuss the basis of 
quasi-static theory. In the latter part, we carry out the molecular dynamics 
simulation of collisions between two identical nanoclusters, and report some 
^ ■ statistical properties of impacts of nanoclusters. 

a 
a ■ 

1 Introduction 

The initial kinetic energy of colliding bodies is distributed into the internal degrees 
of freedom in an inelastic collision. Such a collision is characterized by the restitu- 
tion coefficient e = V /V, where V and V are, respectively, the relative colliding 
speed and the relative rebound speed. Although it is believed that the restitution 
| coefficient e satisfies e < 1 for impacts of macroscopic bodies, the anomalous im- 

pact with e > 1 is possible in some special situations for small bodies. Indeed, the 
prohibition of e > 1 is originated from the second law of thermodynamics [TJ [2], but 
i/-) ■ some terms which disappear in the thermodynamic limit play important roles in 

the description of small systems. It should be noted that the restitution coefficient 
projected into the normal direction of the collision can easily exceed unity in the 
case of oblique collisions. [3] 0] 

The low-speed collisions for macroscopic bodies are believed to be described by 
the quasi-static theory, which is consistent with some experimental results. [5l El 
[7] However, it is not obvious whether the quasi-static theory is applicable to the 
impact of nanoclusters. Indeed, we expect that the effects of cohesive force among 
atoms cannot be ignored for such small systems. Awasthi et al.[8\ reported that 
the dependence of the restitution coefficient e on the impact speed for nanoclusters, 
which contain adhesions, differs from the prediction from the quasi-static theory 
based on their molecular dynamics simulation (MDS). In a recent paper, Brilliantov 
et al. extend the quasi-static theory to the theory of cohesive collisions. [9] 

The physics of nanoclusters is one of hot subjects. The MDS is a standard tool to 
investigate collisions of nanoclusters such as fulleren. Some of such studies focus on 
fragmentations and coalescences after binary collisions of nanoclusters [TUl [TT| PT2] . 
The other studies discuss the collisions of a cluster with a substrate El, the erosion 



process on a diamond surface [T3]. and the fragmentation pattern of clusters [13]. So 
far, we do not know any paper to investigate the effects of thermal fluctuations on 
collisions of nanoclusters except for our preliminary report. [T5] 

In this paper, we perform the MDS of colliding clusters to investigate the ef- 
fect of thermal fluctuations. This is the extension of our previous work. [TH] In the 
first part, we review what model is adequate to describe the collision of nanoclus- 
ters. In the next section, we introduce the generalized Langevin equation and the 
fluctuation-dissipation relation in the first kind. In section 3, we calculate the ve- 
locity autocorrelation function (VACF) of a lattice model to apply it to a system 
of a nanocluster. Then, we justify the quasi-static theory to describe the collisions 
of nanoclusters. In the second part, we show the detailed results of our numerical 
simulations. In section 4, we introduce our numerical model of the MDS. In section 
5, we explain the results of our simulation, which consist of four subsections. In 
the first subsection, we check the relaxation of VACF in our model. In the next 
subsection, we demonstrate sequential snapshots of collisions between two identical 
nanoclusters. In section 5.3, we compare our numerical result of the impact speed 
dependence of the restitution coefficient with the quasi-static theory of cohesive or 
noncohesive collisions. In section 5.4, we show the frequency distribution functions 
of the restitution coefficient and their dependence on cohesive force between atoms. 
We also show the probabilities to appear four categories in our simulation when the 
cohesive parameter is finite. In section 6, we discuss and summarize our results. 

2 The Langevin equation 

It is well known that we can formally rewrite the Newtonian equation of motion 
as the generalized Langevin equation for the 'slow' variable. When we consider the 
motion of colliding a pair of small clusters, it is natural to adopt the relative velocity 
v between the center of mass of each cluster as the 'slow' variable [TE]. We should 
note that the center of mass is characterized by the total mass of one cluster, while 
each element of the cluster can be characterized by the mass for the element. Thus, 
the effective mass of the center of mass is much larger than the mass of the element. 
Thus, the generalized Langevin equation is given by 

where j(t), 0(t), M, and F(t) are the memory kernel, the fluctuating force from the 
fast oscillations, the reduced mass of two clusters, and the systematic force acting 
on the centers of mass, respectively. The fluctuation force 6(t) is believed to be 
unimportant for the impact problem of two clusters. The systematic force F(t) may 
be approximated by the Hertzian contact force. From the fluctuation-dissipation 
relation in the first kind, the Laplace transform j(oo) = J °° dt"y(t)e~ lut satisfies 

(iu + 7M)- 1 = -=■ / dt(v(0) ■ v(t))e-^\ (2) 

1 J —oo 

where T is the temperature, and the Boltzmann constant is set to be unity [17]. We 
should note that 7(0*) can be defined as the usual Fourier transform if we assume 
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7(t) = 7(— t) for t < 0. If the integration of the velocity autocorrelation function 
(VACF), i.e. dt(v(0) ■ v(t)) is finite, j(uj — > 0) is finite. In this case, the 
generalized Langevin equation can be approximated by the Langevin equation with 
the white noise satisfying (9i(t) ■ 6j(t')) = 2T , ySijS(t — t') where 9i(t) is the i-th 
component of 0{t). It is obvious that this Langevin equation is an irreversible 
equation. Thus, the behavior of VACF is the most important to characterize the 
macroscopic dissipation. 



3 The relaxation of the correlation function 

As discussed in the previous section, the relaxation of VACF plays a key role for 
the equilibration process of a system. Let us consider the relaxation of VACF in a 
nanocluster which consists of a regular lattice with equal-mass atoms. When the 
atoms are confined in attractive potential, the excitation from the ground state is 
characterized by the harmonic oscillation in a simple cubic lattice. In this section, 
we demonstrate that VACF of a uniform system in the simple cubic lattice ex- 
hibits the slow relaxation proportional to The analysis presented here is the 
straightforward extension of the one-dimensional cases. [TB] 

Let us consider an infinitely large simple cubic lattice system in which the mass 
points with mass m connecting with the linear spring whose spring constant is k. 
The position of each lattice point can be specified by a set of integer n = (n x , n y ,n z ) 



in this system. Introducing the characteristic angular frequency cuo = yk/m, the 
equation of motion of the deviation from the equilibrium position r n (t) obeys 

6 

r n (t) = -u; 2 (6r n (t) - ^r^*)), (3) 

i=l 

where e x = (1, 0, 0), e 2 = (-1, 0, 0), e 3 = (0, 1, 0), e 4 = (0, -1, 0), e 5 = (0, 0, 1) and 
e 6 = (0,0,-1). 

Let us introduce the lattice Fourier transform and the inverse Fourier transform 

as 

r k (t) = J2 e-^ n r n (t), r n (t) = J dke lk n r k (t), (4) 

where £„ = E^=-oo ££=-00 and J dk = J^dk x J^dk y f^. Thus, the 

equation of motion in the Fourier space is given by 

^(t) = -4c 2 Esin 2 (|)f(t). (5) 

Furthermore, introducing the Laplace transform fk(z) = / °° dte~ zt r k {t), we obtain 

= rf.(0) + *.(0) 

Let us consider the motion of the mass point at the center of the mass in the 
system. From the relation r Q (t) = J T^s^fc^), we reach 

fn(z)- 1 fdk ***(°) + *»(°) (7) 
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where we have used r n (z) = / °° dte zt r n {t). 

Here, VACF at the center of the mass is defined by 

0(i) = (r o (O)-fo(t)>, (8) 

where (• • •) represents the ensemble average over the different initial conditions. We 
assume that the initial condition satisfies 
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(r m (0) • r n (0)> = 0(O)5 m ,» = — $m, n , (r n (0) • ^o(O)) = 0. (9) 

Then, the Laplace transform <j>(z) = J °° e~ zt <f)(t) = z(f (z) ■ ro(0)) of <p{t) satisfies 
~. . _ 1 z(r h {0) ■ f o (0)) + (f fe (0) • f o (0)) 



(2tt) 3 

3Tz 
(2vr) 3 m 



^ 2 + 4^EtiSin 2 (| 

C?fe -. — 

2 2 + 4cu 2 ELsin 2 (| 



(10) 



where we have used (^fe(O) • 'To(O)) = 3T/m. From the inverse Laplace transform of 
<j>(z) we obtain the expression 
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dk 



cos 







2cu t, 









Since the direct integration of ffTTT) is difficult and we are not interested in the 
detailed properties of the lattice, we may introduce the approximation 2f =1 sin 2 ^ ~ 

J2i=i /4 = /c 2 /4 ~ sin 2 |, where A; = \Jj2i=i kf- Once we adopt such an approxi- 
mation, we obtain 
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2TT 2 m Jo 



dkk cos 



2uoot sin • 



k 



(12) 



From the numerical integration of this expression, it is clearly to find cj)(t) ~ l/vt 
(see Fig. CD). Indeed, when we put r = 2co>o^sinfc/2, there is the relation k 2 dk = 
4dk—4:dT /(uJot)+0(k 4 dk). If we ignore the terms of k 4 dk, we obtain the approximate 
relation I(t) = f$ dkk 2 cos[2u;ot sin k/2] ~ 4n J (2u t) — 4sin(2o;ot)/(a;ot). From the 
asymptotic form of the Bessel function, we obtain I(t) ~ J n/ (u t) cos(2uj t — 7r/4) 
for uj Q t ^> 1. This result is essentially the same as that for one-dimensional case 
4>{t) oc Jo(2co>o<0- This time dependence of VACF can be observed in the direct 
simulation of nanoclusters in which each cluster is a 13 layer of the spherical cut 
of face-centered cubic (FCC) lattice , i.e. 682 Lennard- Jones atoms (Fig. 3). The 
setup of our simulation will be explained in the latter part. 

In spite of this extremely slow relaxation under the absence of 7(0), we can 
approximately define the friction constant in the Langevin equation to describe 
the low frequency behaviors. Actually, if we adopt the approximation 0(t) ~ 
a cos(u;ot) / \Ju \t\ with a constant a for large t, substituting this into we ob- 
tain 
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Figure 1: A log-log plot of the decay of 
VACF. Here the time is normalized by 
u Q and the guide line represents 20/ \/t. 
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Figure 2: The comparison between 
VACF of the center of mass of the up- 
per cluster (cross points) and fitting line 
proportional to 1/ \/t, where the time is 
dimensionless time for the MDS. 



for uj <C ujq. Thus, we may approximate 7(t) = /f^ du/2ire lU3t, y(u}) by 

7 (t) ^ n*(*) - ^ (14) 

where 8(t) is Dirac's delta function and Q = T^/uJq/ {\/2am). Thus, the memory 
term can be approximated by 

/' dn 
dt'i{t')v{t') ~ Qv{t) + ^ - S(0)v{t), (15) 
-oo dt 

where the last term can be absorbed in the initial condition. Finally, we obtain the 
effective Langevin equation for the low frequency behavior at t ^ as 

dv n e f 

— = v-\ 1 (16) 

dt 2 2 2M v ; 

which does not have any essential difference from the conventional Langevin equa- 
tion. This may justify to use the quasi-static theory even when we consider a 
collision between nanoclusters of a uniform lattice system. Indeed, once we accept 
to use the Langevin equation, it is straightforward to derive the quasi-static theory 
of macroscopic collisions. [SI El E] 

It should be noted that the motion of the atom at the center of mass can be 
described by the equation of motion for a harmonic oscillator. In order to use eq.(l), 
we need to introduce some tricks, such as the mass difference, the contact with the 
other atoms and the nonlinearity. However, this argument may be instructive to 
understand the basis of the Langevin equation from the mechanical point of view. 



4 Our numerical model 

Let us introduce our numerical model. Our model consists of two identical clusters. 
Each of them is the spherical cut from a 13 layered face-centered cubic (FCC) lattice 
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and consisted of 682 "atoms". When we simulate larger size of nanoclusters, the 
system is fTuidized in the vicinity of surface, while the data for the smaller systems 
strongly depend on the specific orientation of impacts. The details of system size 
dependence of the simulation will be reported elsewhere. 

The clusters have facets because of the small number of "atoms" (Fig. [3]). All 
the "atoms" in each cluster are bounded by the Lennard- Jones potential Ufcj) as 




U{r l3 )=Ae{\-\ -a - , (17) 



where is the distance between two "atoms", i and j. The coupling coefficient 
of the attractive term a is treated as a cohesive parameter between atoms on the 
surfaces of one cluster and those on the surface of another, while the potential act on 
the atoms within the same cluster satisfies a = 1.0. Here, we will consider collisions 
for the control parameters a = and a = 0.2 between different clusters. 0. In 
eq. f|T7|) . e is the energy constant and a is the lattice constant. When we regard the 
"atom" as argon, the values of the constants become e = 1.65 x 10 -21 J and a = 3.4A, 
respectively. |18J Henceforth, we label the upper and the lower clusters as cluster 
C u and cluster C l , respectively. To reduce computational costs, we introduce the 
cut-off length <j c of the Lennard- Jones interaction as a c = 2.5a. 

The procedure of our simulation is as follows. The initial velocities of the "atoms" 
in both C u and C l satisfy the Maxwell-Boltzmann distribution at the initial tem- 
perature T. The initial temperature is set to be T = O.Ole or T = 0.02e in most 
of our simulations. Sample average is taken over different sets of initial velocities 
governed by the Maxwell-Boltzmann velocity distribution for "atoms" . 

To equilibrate the clusters, we adopt the velocity scaling method [T§1 [20] for 
2000 steps in the initial stage of simulations. We have checked the equilibration 
of the total energy in the initial relaxation process. After the equilibration, we 
give translational velocities to C u and C l at the relative separation a c between two 
clusters to make them collide against each other, where the initial colliding speed is 
achieved by the acceleration g = 0.01e/(m<r) from a stationary state. The relative 
speed of impact ranges from V = 0.02-y/e/m to V = 0.07-y/e/m, which are less than 



the thermal velocity for one "atom" defined by JT/m, where m is the mass of the 
"atom". 

Numerical integration of the equation of motion for each atom is carried out by 
the second order symplectic integrator with the time step dt = 1.0 x 10~ 2 a/ Jefm. 
The rate of energy conservation, \E{t) — E$\/\Eq\, is kept within 10~ 5 , where Eq is 
the initial energy of the system and E(t) is the energy at time t. 

We let the angle around z— axis, 8 Z , be 9 Z = when the two clusters are located 
in mirror-symmetric positions with respect to z = 0. In most of our simulation, we 
adopt the data at 9 Z = 0. From our impact simulation for Q\ = 7rz/18 {i = 1, 9) 
at T = 0.02e we have confirmed that the initial orientation does not crucially affect 
the restitution coefficient. 
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5 The results of our numerical simulation 



5.1 The relaxation of velocity autocorrelation function 

At first, we have carried out the contact simulation for two identical nanoclusters 
contacting each other. From our simulation, we verify that the Hertzian contact 
theory can be used without introduction of any fitting parameters [TS]. The details 
will be reported elsewhere. Another purpose of the contact simulation is to check 
whether eq. (fT2l can be used in our system. For this purpose, we make the two 
identical clusters contact each other under the mirror symmetric configuration, and 
equilibrate them at T = 0.03e. After the equilibration, we leave those clusters, and 
record the time evolution of the velocity of 14 atoms near the center of mass of 
the upper cluster, and collect 50 samples with different initial velocities for all the 
atoms, i.e. we average the data under 700 different samples to calculate VACF . 

Figure [2] is the result of VACF near the center of mass of the upper cluster in 
our simulation. The upper envelope line is given by f(x) = 0.0029a; -1 / 2 , which is 
consistent with the theoretical prediction. 

5.2 The collision of two identical clusters 

In this subsection, we show the results of out simulations for colliding two identical 
nanoclusters. We mainly simulate the two cases for the interaction between different 
clusters: the completely repulsive case with a = 0.0, and the weakly cohesive case 
with a = 0.2, where a is the cohesive parameter in eq. (fTTj) between different clusters. 

Let us show the sequential snapshots of two colliding clusters. Figure [3] (a) and 
(b) show the collisional behavior in the case of a = 0.0 and a = 1.0, respectively. It 
should be noted that we demonstrate the case of a = 1.0 to emphasize the difference 
between the noncohesive collision and the cohesive collision. In Fig. [3] (b), we can 
observe the elongation of clusters along the z— axis before the separation, while we 
do not observe any elongation of clusters before the separation in Fig. [3] (a). This 
elongation in Fig. [3] (b) is the result of the cohesive interaction between two clusters. 

5.3 The relations between the restitution coefficient and the 
colliding speed 

Here, we numerically investigate the relation between the restitution coefficient and 
the colliding speed. Figure H] shows the relation between the restitution coefficient 
e and the relative speed of impact V/ \JTJm in purely repulsive collisions with a = 
0. The initial configurations of two colliding clusters are assumed to be mirror 
symmetric. The cross points and error bars in FigJH are, respectively, the average 
and the standard deviation of 100 samples for each colliding speed. From Fig. HJ we 
confirm that the restitution coefficient e decreases with the increase of the colliding 
speed V/^e/m. When the colliding speed is V = 0.02y / e/m at T = 0.02e, the 
average of e becomes 1.04 which is slightly larger than unity. It is interesting that 
our result can be fitted by the quasi-static theory of low-speed impacts 1 — e oc 
V^/ 5 [5J El [7] when the restitution coefficient in the limit V — > is replaced by a 
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Figure 3: Sequential snapshots of two colliding clusters in the cases of (a) a=0.0 and 
(b) a=1.0. 



constant larger than unity. Indeed, the solid and the broken lines in Fig. H] are 

fitting curves of e = a± — 012 (V/ ^Je/rn^J , where a.\ and depend on material 
constants of colliding bodies and T. 

We also briefly discuss the effect of the size dependence on the result. The results 
of our simulation for iV = 433, which is a 11 layered spherical cut of FCC lattice, 
cannot be approximated by the quasi-static theory, where the restitution coefficient 
seems to be almost independent of the colliding speed in the wide range of the 
impact speed. On the other hand, we cannot find any systematic relation between 
the restitution coefficient and the colliding speed in the simulation for iV = 1466 
which is a 17 layered spherical cut of FCC lattice. This can be attributed to the 
melting on the surface of the cluster. The details of the melting properties will be 
reported elsewhere. 

Next, we investigate the weakly cohesive collisions with a = 0.2 between those 
two clusters. Figure [5] shows the relation between restitution coefficient and impact 
speed, where 100 samples are taken for each colliding speed at the initial temperature 
T = 0.02e. When there is the cohesive interaction between two colliding clusters, the 
relation has a peak as suggested by Brilliantov et al. [9]. In the figure, the open circles 
are the numerical results obtained by solving the equation developed by Brilliantov 
et al.[9\. To solve this equation, we evaluate the values 7 ~ 0.026^/em/a from the 
calculation of the attractive interaction between two clusters^] The theoretical result 
in Fig. [5] suggests that the restitution coefficient is insensitive to the colliding speed 



1 The surface tension 7 can be calculated from the attractive potential. The method of our 
evaluation will be reported elsewhere. 
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Figure 4: The relation between col- 
liding speed and restitution coef- 
ficient in the case of a = 0.0 at 
T = O.Ole and T = 0.02e. 
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Figure 5: The relation between col- 
liding speed and restitution coeffi- 
cient in the case of a = 0.2. One 
cluster consists of 682 atoms. The 
cross points and open circles are, 
respectively, the result of our simu- 
lation and the theory in ref.[9]. 



for the large colliding speed, though the restitution coefficient slightly decreases with 
the increment of the colliding speed. 

5.4 The frequency distribution functions of the restitution 
coefficient 

Here, we show our numerical results on the frequency distribution function of the 
restitution coefficient, which strongly depends on the cohesive parameter. Figure 
[6] shows histograms of the restitution coefficients for both the purely repulsive col- 
lisions and the cohesive collisions a = 0.2. When there is no cohesive interaction 
between the two clusters, the frequency distribution function is roughly represented 
by the Gaussian distribution function. On the other hand, the frequency distribu- 
tion function is irregular when the cohesion exists. It is notable that the anomalous 
events for e to exceed the unity becomes rare when there is the attractive interaction 
between clusters, though a few percent of the collisions still exhibit the anomalous 
impacts. This is because two clusters are coalesced with each other in the slow 
impacts. Therefore, the frequency distribution function for a = 0.2 has a steep peak 
near e — 1. 

In Figj6](b), there are the first and the second peaks around e = 0.4 and e = 0.65, 
respectively. The collisional modes observed around these peaks are the rotational 
bounces after the collisions, while the most of bounces are not associated with 
rotations in the vicinity of the third peak around e = 1. It is reasonable that the 
excitation of macroscopic rotation lowers the translational energy to decrease the 
restitution coefficient. 

For cohesive collisions, we can categorize the rebound behaviors of the colliding 
clusters into four patterns (see Fig. [7]): (a) n = (complete adhesion), (b) n > 1, 
(c) n — 1 and e < 1, and (d) n — 1 and e > 1, where n is the number of collisions in 
each impact process. The collision with n > 1 can take place, when the attractive 
interaction between the colliding clusters exists. Indeed, if the rebound speed is not 
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(a) (b) 



Figure 6: Histograms of the restitution coefficients for (a) a = 0.0, V = 0.02ye/m, 

and (b) a = 0.2, V = O.lJe/m. The guide line in the left figure is the Gaussian 
fitting of the data. 



large enough, the rebounded clusters are attracted to have the second collision. We 
call the case with e > 1 and n — 1 the anomalous impact, but there are some other 
characteristic collisions as can be seen in Fig. [71 

Similarly, we categorize the collisions into four groups as a function of the co- 
hesive parameter under the fixing colliding speed V = 0.02-^/e/m (Fig. 8). It is 
obvious that there are two categories, (c) and (d), in noncohesive collisions, while 
the probability to occur (a) or (b) increases as a increases. It is interesting that Fig. 8 
is almost the mirror symmetric one of Fig. 7. This fact suggests that the cohesive 
parameter plays a role of the impact speed. The relation between the impact speed 
and the cohesive parameter will be discussed elsewhere. 
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Figure 7: Probabilities classified by col- 
lision modes observed in cohesive colli- 
sion with a = 0.2. 
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Figure 8: The probability classified by 
four modes at V = 0.02 J e/m. 
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6 Discussion and conclusion 



In this paper, we study collisions of nanoclusters which are thermally activated. We 
also discuss the effects of cohesive force between the colliding two clusters. Although 
the results are preliminary, we believe that our paper includes some potentially im- 
portant results for the nanoscience. Let us briefly discuss our results. An anomalous 
impact with e > 1 occurs with a finite probability even for realistic situations (see 
Fig. 8). This is an important indication, though the cohesive force between the 
colliding clusters suppresses such events in the low speed collisions. We also find 
an interesting similarity in the roles of the impact speed and the cohesive parame- 
ter (Fig. 8). It is more interesting that the cluster is fluidized when the cluster is 
large. There is capillary instability at the surface of the large cluster, because the 
influence of the attractive binding force from the center of mass is weaker, the size 
of the cluster is larger. The quantitative discussion will be discussed elsewhere. 

In conclusion, we study the impact of thermally activated nanoclusters numer- 
ically. We confirm that VACF satisfies <p{t) ~ l/y/i- The restitution coefficient 
seems to be consistent with the quasi-static theory when there is no attractive in- 
teraction between the two colliding clusters, while the restitution coefficient has a 
peak at a finite value of the impact speed, when the attractive interaction exists. 
The anomalous impacts which have e > 1 commonly take place in purely repulsive 
collisions, while such an impacts become rare in cohesive collisions. The frequency 
distribution function satisfies Gaussian for purely repulsive collisions and has some 
peaks in cohesive collisions. 
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